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ABSTRACT 

Context. The scarcity of accretion rate estimates and accretion tracers available for Herbig Ae/Be (HAeBe) stars contrasts with the 
extensive studies for lower mass objects. 

Aims. This work aims to derive accretion rates from the UV Balmer excess for a sample of 38 HAeBe stars. We look for possible 
empirical correlations with the strength of the Ho-, [Oi]6300, and Bry emission lines. 

Methods. Shock modelling within the context of magnetospheric accretion (MA) was applied to each star. We obtained the accretion 
rates from the excess in the Balmer discontinuity, derived from mean values of multi-epoch Johnson's t/6 photometry. The accretion 
rates were related to both mean Ho- luminosities. Ha 10% widths, and [Oi]6300 luminosities from simultaneous spectra, and to Bry 
luminosities from the literature. 

Results. The typical -median- mass accretion rate is 2 x 10"'' Mq yr"' in our sample, 36 % of the stars showing values < 10"' Mq 
yr"', 35% between 10"' and 10"'', and 29% > 10"'' Mq yr"'. The model fails to reproduce the large Balmer excesses shown by the 
four hottest stars (T, > 12000 K). When accretion is related to the stellar masses and luminosities (1 < M./Mq < 6; 2 < L./Lq < 
10'), we derive M.^^^ oc M^ and Lace I-I ". with scatter. Empirical calibrations relating the accretion and the Ha, [Oi]6300, and Bry 
luminosities are provided. The slopes in our expressions are slightly shallower than those for lower mass stars, but the difference is 
within the uncertainties, except for the [O i]6300 line. The Ha 10% width is uncorrelated with Ma^o unlike for the lower mass regime. 
The mean Ha width shows higher values as the projected rotational velocities of HAe stars increase, which agrees with MA. The 
accretion rate variations in the sample are typically lower than 0.5 dex on timescales of days to months. Our data suggest that the 
changes in the Balmer excess are uncorrelated to the simultaneous changes of the line luminosities. 

Conclusions. The Balmer excesses and Ha line widths of HAe stars can be interpreted within the context of MA, which is not the case 
for several HBes. The steep trend relating Ma^c and M, can be explained from the mass-age distribution characterizing HAeBe stars. 
The line luminosities used for low-mass objects are also valid to estimate typical accretion rates for the intermediate-mass regime 
under similar empirical expressions. However, we suggest that several of these calibrations are driven by the stellar luminosity. 

Key words. Stars: pre-main sequence - Accretion, accretion disks - circumstellar matter - protoplanetary disks - Stars: activity - Line: 
formation 

1 . Introduction hot accretion shocks can be modelled, allowing us to explain the 

», ■ . continuum excess and the spectroscopic veiling from tempera- 

^ , Most of the stellar mass acquired during the pre-main sequence ^^^^ differences with the stellar photosphere ( Calvet & Gullbrinl 

■ ■ ■ (PMS) phase comes through accretion from the disk. In addi- ^gg^ r^-^ modelling yields accretion rate estimates, which are 

tion, estimates of the mass accretion rate are needed to anal- ^^^^^ correlate with the strength of spectroscopic features 

yse circumstellar gas dissipation, and therefore to quantify the covering the wavelength rang e from the U V to t he IR (see e.g. 

timescale when planets could be formed m protoplanetary disks iHerczeg & Hillen'M TO iFang etal.l 12001 iRidiaco et all 

^^etospheric accretion (MA hereafter, | Uchida & Shibata| I20III and references therein). Even though the origin of these 

[19851 iKpenigl 1991]; | Shu et al.||1994D is the accepted paradigm correlations is not clear, those spectral lines are used as empiri- 
that explains disk-to-star accretion m classical T-Tauri (CTT) ^^^^^^^^ ^j^^^j^ simplifies the process of estimating accretion 

and lower mass stars. In the MA scenario, the inner disk is j-^^es 
truncated at some point between the stellar surface and the co- 
rotation radius, where matter is accelerated through the mag- Herbig Ae/Be (HAeBe) stars are the massive (1-10Mq) 

netic field lines until it reaches the central stai-. The resulting counterparts of CTTs. How circumstellar matter accretes on 

those objects remains an open issue. For instance, MS stars with 

Send offprint requests to: Ignacio Mendigutia spectral types earlier than ~ A6 (M, > 2 M©) are not expected 

e-mail: Ignacio. Mendigutia@cab . inta-csic . es to have convective sub-photospheric zones generating the neces- 
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sary magnetic fields (ISimon et al.ll2002[ and references therein). 
However, in young stars, convection zones may appear in ear- 
lier spectral types than expected (Finkenzeller 1985). In addi- 
tion, different evidence suggests that MA could be extended to 
the intermediate-mass regime. Although weak in se veral cases, 
magn e tic fields have been detected in some HAeBes (IWade et al.l 
l2007t [Hubrig et alj l2009l) . Spectropolarimetric measurements 
point to MA actiii g in HAe stars ( Vink e t al. 2002, 2003; 
liviottrametaLll2007h . and lGradv et all (1201 Ol) conclude that the 
accretion in these objects goes through high-latitude funnels, 
again suggesting m agnetically controlled accretion (see also 
iBrittain et al.l 12009 ). High-velocity redshifted self-abs orptions 
obser v ed in line profile s of se veral HAeBe stars (Nat ta et aP 
120001: iMoraet all [20021 |2004 point to infalling material at 
close to free-fall velocities. These can be easily explained from 
MA and hardly at all from a competing scenario such as the 
boundary-layer one (see e.g. the review in :Alencaill2007ft . 

MA modelling was caiTied out for a s ample of nine 
intermediate-mass T-Tauri stars bv lCalvet et al.l (l2004t) . showing 
that the correlation between the accretion and Bry luminosities 
can be extended to these objects. The problem o f measuring ac- 
creti on rates for the HAeBe regime was faced bv lMuzeroUe et all 
(12004,. MDCH04 hereafter). They reproduced the line profiles 
of UX Ori by assuming the MA geometry. The absence of sig- 
nificant line veiling and optical excess in most HAeBe objects 
was explained from the similar temperatures characterizing the 
stellar photosphere and the accretion shocks. The excess in the 
Balmer discontinuity was shown to be a valid measurement of 
accretion by MDCH04, providing a calibration between both pa- 
rameters for a typical HAe star. 

The large number of accretion studies for wide samples of 
low-mass stars contrasts with the lack of accretion rate measure- 
ments for HAeBes, which is also reflected in the scarce num- 
ber of empirical tracers of accretion available for these objects. 
Accretion ra te estimates for a wi de sample of HAeBe stars were 
obtained by iGarcia Lopez et al.l t2006) froin the B ry line, us- 
ing the correlation derived by (ICalvet et al.l l2004l) for stars in 
a lower range of stellar temperatures. The recent accretion rate 
estimates by Donehew & Brittain (2011) are based on the cal- 
ibration between the mass accretion rate and the excess in the 
Balmer discontinuity provided in MDCH04, which was mod- 
elled for a specific set of stellar parameters. However, the mea- 
surement of accretion requires modelling observed properties for 
each object individually. The accretion rates can then be com- 
pared with spectral diagnostics to establish possible empirical 
correlations. That comparison could benefit from simultaneous 
measurements, given the variabihty cha racterizing PMS stars 
dEiroa et al.ll2 002'; Mendigu tfa et al.ll201 iL Paper I hereafter). 

Following MDCH04, in this work we apply shock modelling 
within the context of MA, aiming to reproduce the strength of 
mean Balmer excesses derived from multi-epoch photometry of 
a wide sample of HAeBe stars. The accretion rates obtained in 
this way are then related to the mean strength of the Ho- and 
[Oi]6300 emission lines from simultaneous spectra, as well as 
to non-simultaneous Bry luminosities from the literature. In ad- 
dition, the simultaneous multi-epoch data allows us to make a 
first approach on the accretion rate variability and its relation 
with the spectral lines. 

Section |2] summarizes the properties of the sample and the 
data. Section [3] describes the shock model (Sect. 13.11 1 and pro- 
vides the accretion rates estimates (Sect. 13.2b . which are com- 
pared to previous results (Sect. IJXTl i. Section prelates our ac- 
cretion luminosities with those of the Ha, [Oi]6300, and Bry 
lines. The relation between accretion and the Ho- 10% width is 



treated separately in Sect. 14.11 The analysis of the accretion rate 
variability is outlined in Sect. |5] The discussion and conclusions 
are included in Sects. |6]and|2l 

2. Sample properties and data 

Table [T] shows stellar parameters of the 38 stars in the sample. 
Columns 1 to 9 indicate the name of the star, the stellar mass, 
luminosity, effective temperature, radius, surface gravity, age, 
projected rotational velocity (v sin /), and distance. Uncertainties 
can be found in s ome o f the references given in the caption. 
iMontesinos et al.l (l2009h . from where most of the values are 
taken, gives around 12%, 40%, + 150 K, +0.1 dex and 35% 
and for the stellar mass, luminosity, eff'ective temperature, sur- 
face gravity and age, respectively. The low errors, as compared 
to other methods based on spectral types, arise from a detailed 
spectroscopic and photometric modelling, based on the compar- 
ison to Kurucz synthetic spectra and stellar evolutionary tracks. 
The typical uncertainty for vsin/ is 6% (Mora et al. 2001). The 
sample is the same one as was analysed in Paper I. The objects 
are mainly HAe and late-type HBes, as well as ten intermediate- 
mass T-Tauri stars with F-G spec tral types and 1-4 Mq. All the 
stars show IR excess (Merm 2004) indicative of dusty circum- 
stellar environments, and the Ha line in emission, pointing to 
active accretion. 

This work is mainly based on the multi-epoch Ho- and 
[Oi]6300 spec tra from Paper I , and th e multi-epoch UBV pho- 
tometry from Oudmaiier et al. (2001). The spectra were taken 
at mid-resolution (R ~ 5500), and the photospheric contribution 
was subtracted. We selected those spectra and photometry taken 
on the same nights. Therefore, the results in this paper were de- 
rived from (almost) simultaneous spectroscopic and photomet- 
ric measurements, unless otherwise stated. The slit width of the 
spectrograph was 1" projected in the sky, which avoids contam- 
ination from stellar companions in almost all sources. This is 
not the case for the photometry, where the aperture was ~ 14". 
However, most of the stars in our sample with reported mul- 
tiplicity have m uch fainter compan ions in the optical (see e.g. 
the references in lWheelwright et al]|201ft) . We assume that their 
contribution to the emission at short wavelengths is negligible 
for our purposes. 

Columns 10, 11, and 12 of Table[T]list the dereddened mean 
Ha luminosities (<L(Ha)>), the mean Ha widths at 10% of peak 
intensity (<Wi()(Ha)>), and the dereddened mean [Oi]6300 lu- 
minosities (<L([Oi]6300)>). The mean line luminosities from 
Paper I were dereddened using the mean colour excesses in Col. 
13, a total-to-selective extinction ratio Ry = 5 (Hernandez et al. 
2004), and t he extinction law co mpiled in Robitaille et al. 
(i2007if l from iKim et all (11994 and 'inde betouw et al.. (^OOsT 
The mean colour excesses w ere computed using the BV pho- 
tometry in Oudmaiier et al. (2001) and the coiTesponding in- 
trinsic colours in Kenvon & Hai'tmann ( 1995.) . Using a different 
col our to characterize th e reddening in the sample (e.g. E(V-R), 
see ICalvet et~an 12004*) would not significantly affect our re- 
sults. <E(B-V)> is consistent with the colour excess applied by 
IMontesinos et al. 1(2009) to derive most of the stellar parameters 
use d in this work. Regardin g the total-to-selective extinction ra- 
tio, [Hernandez et al.l (l2004l) confirmed that Ry = 5 is more con- 
sistent with the observed properties of HAeBe stars than the typ- 
ical Ry = 3.1 for the interstellar medium. U sing a differe nt Ry 
or extinction law (see e.g. the discussion in ICalvet et al.l 120041) 

' The extinction law can be downloaded from tlie web-based SED 
fitting tool in http://caravan.astro.wisc.edu/protostars/j 
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Table 1. Sample of stars. 



Star 


M. 


L. 


T. 


R. 


g 


Age 


V sin ; 


d 


<L(Hq')> 


<Wio(Ht!')> 


<L([OI]6300)> 


<E(B-V)> 


L(Br7) 




(Mo) 


(Lq) 


(K) 


(R©) [cm s -] 


(Myr) 


(km s" 


') (pc) 


[L©] 


(km s"') 


[L©] 


(mag) 


[Lq] 


HD 31648 


2.0 


21.9 


8250 


2.3 


4.0 


6.7 


102 


146 


-1.42 


595 




0.02 


(?) 


HD 34282 


<2.1" 


5.13" 


9550" 


0.8 


4.9 


> 7.8" 


129 


164" 


-2.82 


487 




0.19 


-4.20' 


HD 34700 


2.4* 


20.0* 


6000* 


4.2 


3.6 


3.4* 


46 


336" 


-2.29 


334 




0.01 


(?) 


HD 58647 


6.0 


911 


10500 


9.1 


3.3 


0.4 


118 


543 


-0.13 


619 


-2.49 


0.13 


-2.08^ 


HD 141569 


2.2" 


22.9" 


9550" 


1.8 


4.3 


6.7" 


258 


99" 


-2.01 


646 


-3.71 


0.09 


-3.99'- 


HD 142666 


2.0" 


17.0" 


7590" 


2.4 


4.0 


5.1" 


72 


145" 


-2.33 


483 


-4.75 


0.26 


-3.53'-3 


HD 144432 


2.0" 


14.8" 


7410" 


2.3 


4.0 


5.3" 


85 


145" 


-1.87 


421 


-4.93 


0.06 


-3.29'-3 


HD 150193 


2.2 


36.1 


8970 


2.5 


4.0 


5.0 


100° 


203 


-1.15 


458 




0.45 


-2.64' 


HD 163296 


2.2 


34.5 


9250 


2.3 


4.1 


5.0 


133 


130 


-1.17 


726 


-4.37 


0.03 


.2.771,2.3 


HD 179218 


2.6 


63.1 


9500 


2.9 


3.9 


3.3 


72^ 


201 


-1.16 


464 


-3.86 


0.08 


-2.743 


HD 190073 


5.1 


471 


9500 


8.0 


3.4 


0.6 


20^^ 


767 


0.06 


378 


-2.49 


0.13 


(?) 


AS 442 


3.5 


207 


11000 


4.0 


3.8 


1.5 


(?) 


826 


-0.15 


646 


-2.42 


0.73 


(?) 


VXCas 


2.3 


30.8 


10000 


1.9 


4.3 


6.4 


179 


619 


-1.43 


672 


-3.48 


0.37 


(?) 


BH Cep 


1.7" 


8.91" 


6460" 


2.4 


3.9 


8.2" 


98 


450" 


-2.34 


705 


-4.25 


0.31 


(?) 


BO Cep 


1.5" 


6.61" 


6610" 


2.0 


4.0 


11.2" 


(?) 


400" 


-2.51 


685 


-3.97 


0.13 


(?) 


SV Cep 


2.4 


37.5 


10250 


1.9 


4.3 


5.2 


206 


596 


-1.33 


731 


-3.20 


0.39 


(?) 


V1686 Cyg 


>3.5" 


257" 


6170" 


14 


2.7 


< 0.2" 


(?) 


980" 


-0.27 


457 


-2.80 


0.63 


-1.773 


RMon 


>5.1" 


2690" 


12020" 


12 


3.0 


< 0.01" 


(?) 


800" 


0.34 


832 


-1.04 


0.70 


(?) 


VY Men 


>5.1" 


15800" 


12020" 


29 


2.5 


< 0.01" 


(?) 


800" 


-0.65 


719 


-0.46 


1.79 


(?) 


51 Oph 


4.2 


312 


10250 


5.6 


3.6 


0.7 


256 


142 


-1.23 


522 




0.03 


-2.68''^ 


KK Oph 


2.2" 


25.7" 


7590" 


2.9 


3.8 


3.9" 


177 


160" 


-2.28 


593 


-3.53 


0.36 


-3.53' 


TOri 


2.4 


50.2 


9750 


2.5 


4.0 


4.0 


175 


472 


-0.88 


680 


-2.95 


0.54 


(?) 


BF Ori 


2.6 


61.6 


8970 


3.3 


3.8 


3.2 


37 


603 


-1.24 


731 


-3.49 


0.15 


-2.92^ 


CO Ori 


>3.6" 


100" 


6310" 


8.4 


3.1 


< 0.1" 


65 


450" 


-0.99 


553 


-2.77 


0.70 


(?) 


HK Ori 


3.0" 


77.6" 


8510" 


4.1 


3.7 


1.0" 


(?) 


460" 


-1.57 


573 


-2.69 


0.37 


-2.92''3 


NVOri 


2.2^" 


21.2^^ 


6750^ 


3.4 


3.7 


A.V 


81 


450' 


-1.97 


583 


-4.81 


0.08 


(?) 


RYOri 


2.5" 


28.2" 


6310" 


4.5 


3.5 


1.8" 


66 


460 


-1.7 


598 


-3.74 


0.49 


(?) 


UXOri 


2.3 


36.8 


8460 


2.8 


3.9 


4.5 


215 


517 


-1.36 


677 


-3.58 


0.17 


-2.80''3 


V346 Ori 


2.5 


61.4 


9750 


2.8 


4.0 


3.5 


(?) 


586 


-1.87 


889 




0.29 


-3.21^ 


V350 Ori 


2.2 


29.3 


8970 


2.2 


4.1 


5.5 


(?) 


735 


-1.39 


724 


-3.26 


0.47 


-2.62^ 


XY Per 


2.8 


85.6 


9750 


3.3 


3.9 


2.5 


217 


347 


-1.12 


728 


-3.29 


0.46 


-2.97^ 


VVSer 


4.0 


336 


13800 


3.2 


4.0 


1.2 


229 


614 


-0.06 


691 


-1.82 


1.04 


-1.34'-3 


CQTau 


1.5'^ 


5" 


6800* 


1.6 


4.2 


7.7* 


105 


130' 


-2.86 


529 


-4.54 


0.25 


-3.96'-3 


RRTau 


5.8 


781 


10000 


9.3 


3.3 


0.4 


225 


2103 


0.04 


681 


-1.47 


0.51 


-1.58' 


RYTau 


1.3^ 


2.30^ 


5770^ 


1.5 


4.2 


6.5° 


55 


134' 


-1.97 


648 


-3.43 


0.37 


(?) 


PX Vul 


1.5" 


5.25" 


6760" 


1.7 


4.2 


14" 


(?) 


420" 


-1.56 


628 


-3.71 


0.45 


-2.80^ 


WW Vul 


2.5 


50.0 


8970 


2.9 


3.9 


3.7 


220 


696 


-0.96 


754 


-3.17 


0.36 


-2.33' 


LkHa 234 


>5.3" 


4680" 


12900" 


14 


2.9 


< 0.01" 


(?) 


1250" 


0.58 


747 


-1.60 


1.02 


(?) 



Notes. Surface gravities and line luminosities are on log scale. means non-detections and "( ?)" unknown values. Un less otherwise stated, 
the stellar masses , luminosities, effective temperatures, gravities, ages, and distances are from iMontesinos et alj i200^) : vsin/ values from 
iMora et al.1 ilOOlh . Ho- and [ Oi]6300 mean values, a nd mean colour excesses from this work. Stellar radius are computed from M. and log 
g when they are available in Montesinos et al.l l l2009h . obtaining equal values if derived from L, and T,. R, comes from these parameters for 
the remaining stars, for which log g is then obtained from M. a nd R,. The photometry f or AS 442 and R Mon was taken from SIMBAD 
ffittp:/ /simbad.u-strasbg.fr/simbad/). Bry luminosities are from ' Garcia Looez et al.' ('2006), 'Brittain et al. ( 2007) and 'Donehew & Brittainj 
j2011l). "Manoi et al. (2006), ^Alons o-Albi e t al. (2009), °Glebocki et al. (200(1), '^ .Guimaraes et al.. C2006) , '^ .Hoffleit & Jascheck C1982.) . 
iMeriiJ i l2004h .' ^iess et al] ( (l999il . lAcke et alj ilOOj) . iBlondel et aP ( l200'6h . 



would not change the major conclusions of this work, given that 
most objects in the sample are not heavily extincted (Ay < 2.5 
magnitudes for 79% of the stars). It is noted that different values 
for Ry wer e used to derive the stella r parameters in Table[T]-e.g. 
Ry = 3.1 in Montesinos et a and a range 1 < Ry < 11 in 

iManoi et"aL ( 2006h -. Section l4l again considers the influence of 
the Ry value on our results. 



Finally, non-simultaneous Bry luminosities from the litera- 
ture (see notes to Table[T]) are given in Col. 14, with mean values 
adopted for the stars where two or three measurements are avail- 
able. For an appropriate comparison, the Bry luminosities were 
rescaled to the same distances used in Paper I to derive the re- 
maining line luminosities (Col. 9). 



3. Accretion rates 

3.1. Description of the model 

HAeBe stars show excess of continuum emission compared 
with MS stars with similar spectral types, in particular at the 
Balmer discontinuity region (see below and iGarrisonI 119781 : 
Donehew & Brittain 2011). The Balmer excesses are modelled 
in this section to provide estimates of the accretion rates. The 
analysis of the origin and strength of the magnetic fields nec- 
essary for driving accretion is beyond the scope of this work. 
We assume the MA geometry, using shock models similar to 
those successfully a p plied f or the lower-mass regime. Following 
iCalvet & GullbringI jT998l) and MDCH04, the total flux per 
wavelength unit emerging from the star is 



F,^fFT' + (l-f)F, 



phot 



(1) 
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where / is the filhng factor that reflects the stellar surface cover- 
age of the accretion columns, F'f'" the flux from the undisturbed 
photosphere, and F"'' the flux from the colum n. The Kurucz 
model corresponding to a given T, and log g (lKuruczlll993h 
is used to represent In turn, the total luminosity of the 

columns is defined from the inward flux of energy carried by the 
accretion columns { T') and the outward stellar radiation below 
the accretion shocks dCalvet & Gullbring|ll998h : 



L™' = F"''A ^{r + FP'"") X A = f Lace + FP'""A 



(2) 



where A - /47rR; is the surface area covered by the shocks. Lace 
= GM,Macc/R» is the accretion luminosity; and ^ - I - R»//?,, 
with Ri the disk truncation radius. The parameter F™' is mod- 
elled as the flux from a blackbody at temperature Ta,i (i.e. F'^"' 
= crT*^i). This is justified for the range of stellar temperatures 
represented here, given that the main contribution to the flux in 
the optical and near-UV comes from the optically thick, heated 
photosphere and that the optically thin pre-shock contribution 
represents less than 1/4 of L™' (see MDCH04). For the model pa- 
rameters used in this work (see below), the blackbody approach 
improves as T* increases, and including the optically-thin con- 
tribution would change our modelled Balmer excesses by less 
than ~ 0.07 magnitudes for the colder stars, which is accurate 
enough for our data. Since the total stellar luminosity can also 
be represented as a blackbody at the stellar temperature, F'''^'" - 
crT^, the second and third terms of Eq.|2]provide 



Finally, from the third and fourth terms of Eq.|2] 

rx/4.Rj.(l-^l^^. 

\ Ri / R, 



(3) 



(4) 



Once T and /?, are fixed, the blackbody temperature char- 
acterizing the flux from the column and the filling factor for 
a given set of stellar and accreting parameters can be deter- 
mined from Eqs. [3]and|4] respectively. We fixed ;F = 10'^ erg 
cm"^ s \ whi ch mostly provides appropriate filling factors (/ 
< 0.1; see e.g. IValenti et al.iri993h and reproduces most of the 
observed Balmer excesses (see Sect. 13.21 1. The use of a differ- 
ent value within the expected interval, 10^'^ < !F(erg cm"^ s ') < 
10'^, does not significantly affect the results (see MDCH04). The 
disk truncation radius s hould be smaller than the co-rotation one 
(Rcor, see Sect. 14.11 and IShu et al ] ll994h . which decreases with 
the stellar rotational velocity: 



^GM.R 



2\l/3 



(5) 



The vsin/ values from Table [T] are used as proxies for v,. 
The disk truncation radius was fixed to 2.5R*, which is the ex- 
pected value for HAeBe stars (MDCH04), or to if this is 
smaller than 2.5Rt. The results do not significantly differ if R,//?, 
is changed up to a factor 2 (Herczeg & Hillenbrand 2008). 

Once Tcoi and / are determined, the expected total flux can 
finally be obtained from Eq. [T] The left-hand panel of Fig. [1] 
shows an example of a photospheric synthetic spectra modified 



^ The use of a family of synthetic models different from Kurucz's 
does not significantly affect our results. The NextGen/PHOENIX mod- 
els ( Hauschildt et al. 1999) provide equal Balmer excesses, within 0.03 
magnitudes, for the typical stellar and model parameters used (see text). 



by the contribution of accretion. This contribution causes the to- 
tal flux to show a filled-in Balmer discontinuity, compared to that 
of the undisturbed photosphere. 

For a given set of stellar and accreting parameters, the mod- 
elled excess in the Balmer discontinuity is defined as AD^ = 
(L'^ - B)pi,oi - (U - B),oiai- This is computed by subtracting the 
U-B colour of the synthetic total flux, which accounts for the 
influence of accretion, from the coiTesponding Kurucz photo- 
spheric colour, which reflects the naked photosphere. The fluxes 
were convolved with the U and B transmission curves pro- 
vided by the Nordic Optical Telescope (left panel of Fig. [TJ, 
since the photometric observations were carried out with the in- 
strument TurPol on that telescope dOudmaiier et al ]|200D. That 
the photometric points are slightly displaced from the spec- 
tra is caused by the strong gradients in the fluxes along the 
transmission bands, specially in U. The difference between the 
Kurucz colours computed in this way and the intrinsic colours in 
iKenvon & HartmannI (Il995h is < 0.05 magnitudes for the stellar 
parameters considered here, being typically around 0.01 magni- 
tudes. This accuracy is enough for our purposes and means that 
the modelled excesses can later be compared to the observed 
ones, described below. The right-hand panel of Fig. [T] shows 
a representative subset of the AD^-Macc calibrations we have 
modelled. The curves rise with the corresponding filling factors, 
following Eq. |4] and are limited by / < 1 . There is a strong de- 
pendence on the stellar temperature and surface gravity that must 
be considered when a mass accretion rate is associated to a given 
ADb value. 

From the observational perspective, the mean excess in the 
Balmer discontinuity for a given star is <ADb> - {U - B)o 
- <U - B>dered, wh e rc (U - B)o is the intrinsic colour from 
iKenvon & HartmannI (Il995h . and <U - B>dere d the dereddened 
mean colour, obtained from the observations in lOudmaijer et al.] 
(2001). The dereddening was applied using the <E(B-V)> val- 
ues in Table [1] Ry - 5, and the same extinction law mentioned 
in Sect. |2] The top panel in Fig. |2] shows the mean observed 
and intrinsic B - V colours from which the dereddening was 
applied. The hottest stars in our sample are the most reddened. 
The bottom panel shows the <U - B>dered and (L'^ - B)o colours 
from which the <ADb> values are obtained. The hottest stars are 
also the ones showing the largest Balmer excesses. It worth men- 
tioning that a slight variation in the adopted value for the stellar 
temperature translates into a strong change in <ADb>, for T, 
< 6400 K. This is caused by the steep dependence of both the 
U-B and B - V intrinsic colours on that range of stellar tem- 
peratures (e. g. differences up to 0.2 mag nitudes for variations of 
+ 150 K, see IKenvon & HartmannI 19951) . In turn, this could sig- 
nificantly change the accretion rate estimate (e.g. a variation of 
0.2 magnitudes in the Balmer excess translates into a change up 
to two orders of magnitude in Mace for T* = 6500 K; see right 
panel of Fig. (T). The dependence of the intrinsic colours on T, is 
smoother for almost all the stars studied here (Fig.|2]l. The typi- 
cal uncertainty for <ADb> is between 0.05 and 0.1 magnitudes 
in our sample. This was estimated by considering both uncer- 
tainties in the photospheric colours, from variations in T, of + 
150 K, and the photometric uncertainties, which are the dom- 
inant contribution for most stars. A typical uncertainty of 0.07 
magnitudes is adopted in this paper for <ADb>, but that strictly 
depends on each star and could increase for the coldest objects 
in the sample. 

Individual calibrations, similar to the examples in the right- 
hand panel of Fig. [T] were constructed for each star from the 
corresponding M,, T*, R», and log g parameters in Table [1] In 
this way, a mass accretion rate was assigned to the modelled 
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Fig. 1. (Left): Photospheric flux (grey line), contribution from accretion (fF" , dashed line) and total flux (solid black line) for M» 
= 2.5 M0, = 2.6 R0, T. = 9000 K, ^ = 10'^ erg cmr'^ s"' {T,oi = 12470 K) and Rj = 2.5R,. The mass accretion rate is 5 x 10"^ 
M0 yr"\ (/ = 0.084). All fluxes are computed at the stellar surface. The transmission curves of the U and B filters are plotted in 
violet and blue, respectively. The convolution of the fluxes with the filter responses provide the synthetic photometric points used 
to obtain the modelled photospheric and total colours (from the shaded and filled linked circles, respectively). (Right): Predicted 
excess in the Balmer discontinuity as a function of mass accretion rate for T» = 6500, 9000, and 12500 K (red, green, and blue lines) 
and log g = 4.0 and 3.0 (solid and dashed lines). T - 10'^ erg cm"^ s ' and - 2.5Rt for all cases. 



Balmer excess that matches the observed one. Then L^cc was 
obtained from Mace, M,, and R,. This method does not need dis- 
tance calibrations since it is based on colour excesses. In ad- 
dition, the stellar parameters of most stars in the sample were 
obtai ned from a distance-independent method dMontesinos et alJ 
12009 ). However, distance uncertainties are indirectly introduced 
through Rt -and log g- when this was derive d from a value of L, 
based on a distance (e.g. lManoj et al.ll2b06l) . 

3.2. Results 

Table |2] shows the observed (mean) excesses in the Balmer 
discontinuity and the accretion rates and best model parame- 
ters that reproduce them. The uncertainties for M-^^c were es- 
timated by varying AD^ by ±0.07 magnitudes in the models, 
which is the typical uncertainty for the observed excess (Sect. 
13. 11 1. Upper limits are provided when <ADb> is lower than that 
value. The uncertainties for Lace also consider the typical error 
in the M,/R, ratios. In addition to the sources of error discussed 
in these sections, we refer th e reader to the broad analysis in 
iHerczeg & Hillenbrand! (l2008h on how different uncertainties of 
the stellar and model assumptions could affect the accretion rates 
estimates. 

The high <ADb> value shown by VV Ser cannot be repro- 
duced using the typical value - lO'^ erg cm"^ s"'. Increasing 
;F by one order of magnitude provides Mace ~ 10""* Mq yr"', 
but the filling factor is almost 100%. Further increases in f pro- 
duce lower filling factors, but the corresponding accretion rates 
also increase. Another possibility is to increase the disk trunca- 
tion radius to a value significantly higher than the co-rotation 
one (~ 1.7 R, for VV Ser). For = 5.5 R. and T = lO'^ erg 
cm"^ s the Balmer excess of VV Ser can be reproduced from 
Mace ~ 5x10"^ Mq yr"' (f = 0.99). Similarly, we find it difli- 



cult to reproduce the strong Balmer excesses shown by R Mon, 
VY Mon, and LkHa 234 from the model used in this work. 
For the highest accretion rates possible (~ 10^^-10"' Mq yr"', 
with » 10'^ erg cm"^ s ' and / ~ 1), their modelled spec- 
tra are featureless and provide Balmer excesses that are ~ 50% 
lower than observed. Again, a combination of large disk trun- 
cation radius, well above the typical 2.5 R, for HAeBe stars, 
'F » 10^^ erg cm"^ s and f ~ 1 could reproduce the Balmer 
excess of these objects from lower accretion rates (~ 10"^ Mq 
yr~'). However, the use of those values would be arbitrary and 
not supported by theory or observations (see MDCH04). In par- 
ticular, a disk truncation radius larger than the co-rota tion one is 
in principle not possible in the MA scenario (see e.g. I Shu et al] 
1994). Temperatures and extinction properties that are different 
than the ones used in this work would make their Balmer ex- 
cesses in Table|2]wrong, and therefore useless in our fit attempts. 
However, It is noted that the four problematic objects not only 
show the highest <ADb> values in our sample, but are also the 
hottest (T, > 12000K, i.e. spectral-types earlier than ~ B8) and 
most heavily extincted sources. There is consensus in the liter- 
ature that these objects are early-type HAeB e stars (see e.g. th e 
compilation of spectral types in Table 6 of iMora et alJ l200"ll) . 
Our previous results, indicating that the variability of the Ha 
line is clearly lower for the massive early-type HAeBe stars, al- 
ready pointed to different physical processes operating in their 
circumstellar environments (see Paper 1). The lack of success in 
reproducing the Balmer excesses could be pointing to a change 
in the accretion paradigm for the HBe regime (see also Sect. 
14.11 ). No result is included in Table |2] for the four objects dis- 
cussed, which would require specific analysis. The calibrations 
with the line luminosities included in Sect. |4] and the stellar pa- 
rameters used in this work provide mass accretion rates (in Mq 
yr-i) of - 10-5-6x10-^ for VV Ser, 3x10 ^-3x10-5 for R Mon, 
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Table 2. Observed mean excess in the Balmer discontinuity, ac- 
cretion rates, and model parameters. 



6000 7000 8000 9000 10000 11000 12000 13000 14000 
T. (K) 

Fig. 2. Observed mean colours in B - V (top) and U -B (bottom) 
as a function of the stellar temperature. The colours in the bot- 
tom panel are deredenned using the excesses derived from the 
top panel. T he dashed lines indica t e the c orresponding intrinsic 
colours from lKenvon & Hartmannl(ll995h . The typical eiTorbars 
considering the photometric uncertainties are shown. 



3x10-^-7x10-*' for VY Mon, and exlO^^-SxlO^^ for LkHa 
234. These estimates must be viewed with caution, given the 
discussion above and their being based on relations derived for 
colder stars. 

The typical (median) accretion rate for the remaining objects 
is 2 X 10"^ M0 yr"', 36% of the stars showing Mace < 10"^ 
Mq yr"', 35% between 10"^ and lO"*" Mg yr"', and 29% larger 
than 10"* M0 yr"'. This distribution is practically the same as 
predicted from Fig. 8 of MDCH04. 

Figure [3] (left panel) shows the mass accretion rates against 
the stellar masses. The stars in our sample follow a steeper trend 
than the lower mass o bjects , for which Mace M^; (see e.g. 



MuzeroUe et al ] I2003L I2005L iMohantv et aP l2005t G^ta et alJ 
20061) . We find Mace °^ M^-^ for our sample (Pearson's corre- 



lation coefficient 0.74). The exponent decreases to 4.6 if the up- 
per limits for the accretion rates are discarded (Pearson's cor- 
relation coefficient 0.91). Our data is therefore reasonably well 
fitted to Mace °^ M^. This trend is driven by the most massive 
sources and is most probably caused by t he mass-age distribu- 
tion characterizing HAeBe stars (see e.g. van B oekel et al.l2005t 
iGarcia Lopez et ai.i,2006) . This distribution is shown in Fig. |4] 
for our sample. The stars with M* > 2.5M0 are younger than 
~ 2 Myr, and the objects with lower stellar masses are older. 
Following van Boekel et al. (2005), the lack of old HAeBes 
around 2.5-6 stars is likely caused by the ir faster evolution to 
the M S, which is supported by observations dRoccataghata et all 
l20Tll) . The scarcity of young HAeBes around l-2.5Mo stars is 
most probably explained by the fact that they become optically 
visible later in their evolution. The massive young stars in our 
sample show the strongest accretion rates, as expected. The mass 
accretion rates of the less massive older HAeBe stars are orders 



Star 


<ADs> 


log Mace 


log Lace 


Ri 


T cot 


/ 




(mag) 


[M© yr"'] 


[Lq] 


(R.) 


(K) 


(%) 


HD 31648 


0.05 


<-7.23 


<0.20 


2.5 


12215 


1.1 


HD 34282 


0.06 


<-8.30 


<-0.40 


2.5 


12695 


2.2 


HD 34700 


0.00 


<-8.30 


<-1.05 


2.5 


11730 


0.02 


HD 58647 


0.18 


-4.84+0.22 


2.47 +0.23 


2.1 


13140 


12 


HD 141569 


0.09 


-6.89±0.40 


0.70 +0.40 


1.5 


12695 


3.6 


HD 142666 


0.18 


-6.73±0.26 


0.69 +0.27 


2.5 


12030 


3.2 


HD 144432 


0.06 


<-7.22 


<0.21 


2.5 


11990 


1.1 


HD 150193 


0.29 


-6.12+0.14 


1.33 +0.15 


2.5 


12460 


13 


HD 163296 


0.02 


<-7.52 


<-0.03 


2.2 


12570 


0.61 


HD 179218 


0.02 


<-7.30 


<0.14 


2.5 


12670 


0.60 


HD 190073 


0.22 


-5.00±0.25 


2.29 +0.26 


2.5 


12670 


12 


AS 442 


0.48 


-5.08±0.11 


2.37 +0.12 


2.5 


13405 


56 


VX Cas 


0.22 


-6.44±0.22 


1.16 +0.23 


2.0 


12895 


13 


BH Cep 


0.01 


<-8.30 


<-0.94 


2.4 


11800 


0.07 


BO Cep 


0.21 


-6.93+0.28 


0.45 +0.29 


2.5 


11825 


2.8 


SVCep 


0.24 


-6.30+0.20 


1.30 +0.21 


1.8 


13015 


14 


VI 686 Cyg 


0.12 


-5.23+0.41 


1.66 +0.41 


2.5 


11755 


0.87 


RMon 


0.76 


(?) 


(?) 








VY Mon 


1.22 


(?) 


(?) 








51 Oph 


0.00 


<-7.00 


<0.37 


1.3 


13015 


0.11 


KKOph 


0.36 


-5.85+0.15 


1.51 +0.16 


1.6 


12030 


9.4 


T Ori 


0.09 


-6.58+0.40 


0.90 +0.40 


1.8 


12780 


3.7 


BFOri 


0.00 


<-8.00 


<-0.60 


2.5 


12460 


0.09 


CO Ori 


0.38 


-5.20+0.18 


1.93 +0.19 


2.5 


11775 


4.5 


HK Ori 


0.57 


-5.24+0.12 


2.13 +0.13 


2.5 


12300 


31 


NV Ori 


0.00 


<-8.30 


<-0.99 


2.5 


11850 


0.03 


RYOri 


0.11 


-6.65+0.33 


0.59 +0.33 


2.5 


11775 


0.75 


UX Ori 


0.05 


<-6.77 


<0.63 


1.5 


12285 


1.1 


V346 Ori 


0.25 


-5.99+0.17 


1.47 +0.18 


2.5 


12780 


15 


V350 Ori 


0.20 


-6.66+0.24 


0.82 +0.25 


2.5 


12460 


5.0 


XY Per 


0.16 


-5.86+0.20 


1.57 +0.21 


1.5 


12780 


7.5 


VV Ser 


0.54 


(?) 


(?) 


1.7 






CQTau 


0.02 


<-8.30 


<-0.84 


2.5 


11860 


0.21 


RRTau 


0.37 


-4.11+0.16 


3.18 +0.17 


1.3 


12895 


27 


RYTau 


0.24 


-7.65+0.17 


-0.22+0.18 


2.5 


11700 


0.98 


FX Vul 


0.33 


-6.72+0.16 


0.73 +0.17 


2.5 


11850 


7.2 


WW Vul 


0.08 


-6.38+0.70 


1.05 +0.70 


1.5 


12460 


2.7 


LkHa 234 


0.88 


(?) 


(?) 









Notes. A typical uncertainty of 0.07 magnitudes is adopted for all 
<ADg> values. The symbol (?) refers to the four objects for which 
strong difficulties are found to reproduce their Balmer excesses from 
our model. 



of magnitude lower. These facts would make ffie Macc-M» trend 
very steep. 

The right-hand panel in Fig. [3] shows the accretion lumi- 
nosities against the stellar ones. The stars in our sample have 
0.01L» < Lace ^ L f Both hmits are similar t o those reporte d 
for lower mass stars jClarke & Pringlell2006-.Timng et al.l2008l) . 
The lower limit is the observational detection ffireshold for Laee- 
However, although the trend between Laee and L, shown by low- 
mass stars (Laee °^ ^\'^\ scc e.g. Tilling et al. 2008) is roughly 
followed by our sample, our best fit (Lace °^ L!^) could indicate 
a decrease in the exponent for the intermediate-mass regime. In 
fact, a direct visual inspection of the right-hand panel in Fig. [3] 
reveals that there is an apparent decline in the slope of the lower 
envelope of the Lacc-L» trend. That decrease could start at lower 
stellar luminosities than those analysed in this work, around 
= Lo. 
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Fig. 3. Mass accretion rate vs stellar mass (left) and accretion luminosity vs stellar luminosity (right). Our data are indicated with 
circles, and triangles for the upper limits on the accretion rates. Crosses are data from the literature (those with vertical bars are 
upper limits), including low-mass stars from different star-forming regions (see Natta et al. 2006, and references therein). Dotted 
lines indicate Mace + 1 dex (left) and Lacc/L* = 0.01, 0.1, 1 (right). The dashed line on the right side is the best fit for low-mass 
stars (Lace L'-^), and the solid lines are the best fits for our sample (Mace °c and Lace L'^). 




Fig. 4. Stellar age versus stellar mass for the stars in the sample. 
Error bars are plotted when provided in the references of Table[T] 
The dashed line indicates the theoretical time to reach the main 
sequence (iTavler..l994) . 



3.2.1. Comparison with previous results 

iGarcia Lopez et ail (l2006h obtained accretion luminosities for 
a sample of HAeBe stars from the empirical calibration with 
the Biy luminosity for intermediate- mass T Tauri stars (log 
(Lacc/L©) = 0.9xlog (LB,y/Lo) + 2.9: ICalvet et"alll2004 . The 
non-photospheric Bry equivalent widths were transformed to lu- 



minosities from the K magnitudes and distances in their Table 1 . 
The mass accretion rates were then derived from their values for 
Mt and R,, using the same L^cc-M-^cc expression as in this work 
(Sect. l3.lT ). Figure[5]shows our mass accretion rates against those 
from Garcia Lopez et al. (2006) for the objects in common. The 
difference between both determinations decreases once the Bry 
luminosities are reseated to our distances (see Table [TJ and the 
coiTespoding accretion luminosities are transformed into mass 
accretion rates from our values for M* and R,. The grey symbols 
in Fig.|5]are ob tained in this way, repre senting the estimates from 
the method in iGarcia Lopez et alJ (120061) reseated to our values 
for d, Mt, and R,. The new estimates show a linear correlation 
with our results (Pearson's coefficient 0.65) with slope LO + 0.3, 
and displaced in the ordinate axis by a factor ~ 0.4. This agrees 
with the analysis in the following section. It must be noted that 
the ste llar temperature for several objects in IGarcia Lopez et al.l 
(I2006h falls well above the upper limit from which the Bry rela- 
tio n inlCalvet et alJ (12004 was obtained -around 6100 K-. 

iDonehew & Brittaiiil (1201 1ft have derived mass accretion 
rates for a sample of HAe Be stars from the calibration with the 
Balmer excess provided bv lMuzerolle et all (l2004l) . The Balmer 
excesses were obtained fitting both Kurucz and standard main- 
sequence spectra to UV-optical spec tra. Figure[5]shows our mass 
accretion rates versus the ones from lDonehew & Biittainl (1201 ih 
for the stars included in both samples. Both estimates are equal 
within ± 1 dex for almost all objects in common. Variability 
could account for differences up to ~ 0.5 dex (see Sect. |5]l, but 
the scatter can be explained also from other factors. First, the 
method used to derive the Balmer excess is different in both pa- 
pers. Second, the stellar parameters differ significantly for sev- 
eral stars (e.g. A T, > 10000 K for VI 686 Cyg). Finally, the cal- 
ibration used bv lDoneh ew & Bri ttainI (1201 1 ) was in many cases 
applied to stars with stellar parameters significantly different 
than those from which the calib ration was made (spec tral type 
A2, M. = 3 Mq and R. = 3 Rq iMuzeroUe et alJl2004 . As dis- 
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log M [M yr ], literature 



Fig. 5. Mass accretion rates der iv ed in this work against those 
derived by Garcia Lopez et alJ (l2006l) (filled symbols) and 
iDonehew & BrittainI (1201 1 ) (open symbols). The dashed line in- 
dicates equal values and the dotted lines +1 dex. The typical 
minimum uncertainty for the estimates in iGarcia Lopez et alJ 
(l2Q06h is + 0.5 dex. Triangles are upper li mits. Grey symbols 
repres ent the estimates using the method in iGarcia Lopez et all 
from our distances and stellar parameters, and the solid 
line is the best fit. 



rameters in ICalvet et al.l (l2004l) . The different extinction treat- 
ment is also relevant, but other factors such as the use of differ- 
ent values for T and have less influence. The mass a c cretio n 
rate for the other object in common with ICalvet et al.l (l2004l) . 
RY Tau, is a factor ~ 3 higher than our estimate, but the values 
overlap (~ 10"** Mq yr"') if considering the uncertainties from 
both papers. Apart from the practically equal values for U - B 
and B - V, again within 0.05 magnitudes, the stellar parameters 
adopted in both works are more similar than for CO Ori. 

Given the relevance of the stellar characterization to derive 
the accretion rates, we notice that most of the stellar parameters 
in Table [T] were carefully obtained by Montesi nos et al.l (l2009h 
in a paper specifically devoted to that task, which takes much of 
the photometry used in this work into account. 

4. Accretion tracers 

Figure |6] shows our accretion luminosities in Table |2] against the 
Ha, [Oi]6300 and Bry luminosities in Table [T] The best least- 
squares fits provide the following empirical calibrations for our 
sample (1 < MJMq < 6; 6000 < T,(K) < 1 1000): 



(6) 



log 7^ = 2.28(+0.25)-i- 1.09(+0.16)xlog --^ 



log ( ^ I = 4.80(+0.50) + 1.13(+0.14) x log 



log 



Lo 



3.55(+0.80) + 0.91(+0.27) x log 



-[OIJ6300 I 

L© / 



L^ 



(7) 



(8) 



These empirical expressions are very similar to those pro- 
vided for low-mass CTT stars (Muz erolle et al.l Il998at iDahml 
2008; Herczeg & Hillenbrand 2008): 



log 



2.27(+0.70) + 1.31(+0.16) x log 



-0 



Lac 



log(^) = 6.50(+2.18)+ 1.67(+0.28)xlog 



-[OIJ6300 \ 



(9) 



(10) 



cussed in Sect. 13.11 the ADg-Macc calibration differs strongly 
depending on the stellar parameters considered (right panel in 
Fig.D. 

Th e mass accretion rate provided for CO Ori in ICalvet et al.l 
(12004') was derived from a similar method to the one used in this 
paper, but it was a factor ~ 70 lower. This is explained mainly 
from differences in the stellar parameters adopted. Both L, and 
M» are higher in the present work by a factor ~ 45 and 1.4, 
respectively. From L, and T, (see below), our R, is a factor 2 
larger. The values adopted for R, and M» affect those for the 
filling factor and the blackbody temperature representing the ac- 
cretion shock. The values for the surface gravities and effective 
temperatures are log g = 3.6 against our 3.2, and T, = 6030 
K against our 6310 K. This difference of almost 300 K in that 
range of stellar temperatures, along with the slight differences in 
the U -B and B-V colours -around 0.05 magnitudes-, translates 
into a significa nt change in the B almer excess -0.24 magnitudes 
from the data in lCalvet et al.l (12004 ) against our 0.38 magnitudes; 
see the discussion in Sect. 13.1k The higher mass accretion rate 
we provide agrees with the one expected for a star that here is 
considered hotter, with a lower value for log g and with more 
Balmer excess (see right panel in Fig. [TJ, but our model yields 



log I — U 4.43(+0.79) + 1 .26(+0.19) x log — 



(11) 



There is an apparent decrease in the slope of the calibrations 
for the HAeBe regime, but firm conclusions can only be stated 
for the [O i]6300 line; the slope in Eq.|7]is shallower than that in 
Eq.[TO]above the uncertainties. However, the Ho- and Bry slopes 
for inter mediate-mass T-Ta uri stars are 1.18+0.26 (Dah_m 20081) 
and 0.9 (ICalvet et al.|| 2004]). respectively, i .e. closer to our esti- 
mates (see also Donehew & Brittainll201 lb . This could be indi- 
cating that the slopes of the empirical calibrations decrease for 
higher stellar masses. The scatter for a given line luminosity is 
comparable for classical T-Tauri and HAeBe stars. The slightly 
higher dispersion for the Bry line could likely be diminished us- 
ing simultaneous measurements for the accretion and line lumi- 
nosities. 

As introduced in Sect. |2] a different value for the total-to- 
selective extinction ratio would influence only the Balmer ex- 
cesses of the most heavily extincted sources in our sample. These 
tend to be the early-type HAeBe stars, which are the strongest 
accretors (see Fig. |2] and Sect. 13.2b . Using Ry < 5 would make 
their accretion rate estimates lower. Therefore, the Lacc-L* and 



the same Mace ~ 10 ^-10 ** Mq yr ' for CO Ori using the pa- Lacc-L/me trends for our sample in the right-hand panel of Fig. |3] 
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Fig. 6. Accretion luminosity against the Ha (left), [Oi]6300 (middle) and Bry (right) luminosities . Our results are represented 
with circles and triangles for the upper limits on the accretion rates. Our accretion and Bry luminosities were not derived from 
simultaneous measurements. Error bars for the line luminosities indicate the uncertainties for the non-variable objects and variability 
amplitudes for the variable ones. Our best fits are pl otted with solid lin es. Data for lower mass stars (crosses) were taken from 
iFang et al.l (l2009l) : lHerczep & Hil lenbrand ( 200l), and lCalvet et al.1 (|2004 . for the left, middle, and right panels, respectively. Upper 
limits are indicated with bars, when provided. Dotted and dashed lines are the empirical calibrations for Ha and Bry in low and 
intermediate-mass CTT, respectively (see text). 



and in Fig. |6]would flatten slightly. This would enhance the dif- 
ference between the relations for the low-mass and the HAeBe 
regimes. In addition. The M^cc-M, relation in the left-hand panel 
of Fig. [3] would show a flatter slope, but the trend would be still 
steeper than for low-mass stars, even when using Ry - 3.1. 

4.1. Accretion and the Ha 10% width 

The Ha width at 10% of peak intensity is widely used as an em- 
pirical accretion t racer for low-mass stars ([White & Basrill2003t 
iNatta et al.ll2004t lJavaward hana et"ani2006 l). FigurelTlshows the 
mass accretion rates agains t the Ha 10 % widths. The trend fol- 
lowed by low-mass objects dNatta et a l. 2004) breaks for HAeBe 
stars, whose <Wio(Ha)> values show no correlation with log 
Mace- This was expected since the mean values and the relative 
variabilities of the Ha equivalent and 10% widths are not corre- 
lated in our sample (see Paper I). We conclude that <Wio(Ha)> 
can usually not be used as a direct tracer of accretion for HAeBe 
stars, unlike in the lower mass regime. We argue in the follow- 
ing that one main reason for this discrepanc y is that the typ- 
ically hig h rotation rates of massive stars ( Fi nkenzelled 1 1 985t 
iMora et al. 2001) influence their <Wi()(Ha)> values. We sug- 
gest that this influence can be interpreted in the context of the 
MA scenario. 

Figure [8] rela tes <Wio(Ha)> and v sin / for the stars in the 
sample (see also iJavawardhana et alJl2006l) . Although there is 
a scatter at low velocities, all the rapid HAe rotators (vsin/ > 
150 km s"') have large Ha widths (<Wi()> > 600 km s"'). Our 
spectroscopic measurements refer to the non-photospheric con- 
tribution (see Paper I). Therefore the maximum Ha width must 
be influenced by the circumstellar gas rotating at the maximum 
speed, which, for Keplerian rotation, is located in the inner disk. 
The trend in Fig. [8] suggests that the inner disk gas traced by Ha 
is coupled to the rotating stellar photosphere. Magnetospheric 
channels could be the origin of this link. 



The influence of the stellar rotation on the Ha line broad- 
ening can be qualitatively modelled from MA. The MA model 
described in Muzerolle et al. (2001) was successfully used to re- 
produce several spe ctral lines of accreting stars with differen t 
masses (MDCH04: fCalvet et all [20041: iMuzerolle et alJ l2005h . 
That assumes a simple dipole geometry, where the Ha line is 
generated from free-falling gas that follows the magnetic lines 
connecting the magnetosphere in t he inner disk an d the stellar 
surface. We applied the model from lMuzeroUe et"an (2001) with 
input parameters of a typical HAe star Apart from other effects 
described in that paper, the Ha line broadens whenever the stel- 
lar rotation is increased. The strength of the line broadening also 
depends on additional parameters such as the size of the mag- 
netosphere, the gas temperature, the mass accretion rate, and the 
inclination to the line of sight. The top panel of Fig. |9] shows 
the synthetic profile best fitting the mean Ha line shown by BF 
Ori. This fit adds to that for UX Ori in MDCH04, again exem- 
plifying that MA modelling is able to reproduce the line profiles 
of HAe stars with reasonable accuracy. The modelled profile for 
BF Ori includes the corresponding rotation rate (v sin / = 37 km 
s^^^Mora et al. 2001), and is broadened with respect to the same 
synthetic profile that does not include rotation. The bottom panel 
of Fig. |9] plots some other examples of MA synthetic Ha pro- 
files showing line broadening as the stellar rotation increases. 
All profiles were modelled from the same parameters as in the 
top panel, a magnetosphere size of 2.8-3 R» and different incli- 
nations to the line of sight, as indicated (closer to pole-on for the 
strongest lines and closer to edge-on for the weakest). 



In addition, the trend in Fig. [8] is quantitatively consistent 
with MA operating in HAe stars. When assuming this scenario, 
the region where the gas is channeled through the magnetic ac- 
cretion flows has to be locate d between the st ellar photosphere 
and the co-rotation radius (e.g. lShu et alJl994l) . At this point the 
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Fig. 7. Mass accretion rate against the Ha 10% width. Our re- 
sults are represented by circles, and triangles for the upper limits 
on the accretion rates. Error bars for the line widths indicate the 
uncertainties for the non-variable objects and variability ampli- 
tudes for the variable ones. Data for lowe r mass stars are also 
included (crosses; from the compilation of Natta et al. 2004), as 
well as the empirical expression relating both parameters in that 
stell ar regime; log M arc ~ -12.9 + 9.7x 10"^Wio(Ha) (dotted 
line; lNatta et"al]|2004l) . 



stellar and the disk angular velocities (assuming Keplerian rota- 
tion) are equal, deriving: 



(12) 
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Fig. 8. Ha mean line width at 10% of peak intensity (from Paper 
I) vs the projected rotational velocities (circles). Lower limits for 
<Wi()(Hq')> are marked with vertical bars. Open circles indicate 
HBe stars -51 Oph the object with the highest value of vsin/-. 
Filled circles represent the stars in our sample with later spectral 
types. The theoretical lower limit estimates for each object are 
derived from Eq.[T4] and are plotted as open triangles (see text). 



with Vgas the azimuthal velocity of the gas at the co-rotation ra- 
dius. Since v» > vsin/, and assuming that the non-photospheric 
Ha line width reflects the velocity of the gas at some point be- 
tween R» and Rc 
obtain 



in the sense that Vgas < <Wio(Ha)>/2, we 



< Wio(Ha) > R. 
2v sin / 



(13) 



For the <Wio(Ha)> and vsin / values plotted in Fig.[8l Rcor 
of HAeBe stars is typically in the range ~ 1-5 R*, which agrees 
with the values expected for these objects (see Eq. |5]). We note 
that Rcor < 2 R, for all objects with v sin / > 150 km s Simple 
dipole geometries would need modifications for such small mag- 
netospheric sizes, in order to reproduce the large line strengths 
of HAe stars from MA (see MDCH04). 

Using Eqs. Island [T3l a lower limit for <Wio(Ha)> can be 
estimated as a function of the stellar parameters: 



< Wio(Ha) >„ 



GM» vsin/ 



1/3 



(14) 



We note that this expression is derived by assuming that the line 
broadening is determined by Keplerian gas located below the 



co-rotation radius, which is a major requirement in the MA sce- 
nario. 

Figure|8]overplots our <Wio(Ha)>„„„ estimates from Eq.[T4l 
They reproduce the overall trend followed by our <Wio(Ha)> 
measurements. The largest differences between the observed and 
the estimated values occur for the objects with the lowest values 
of vsin/. This can be partly attributed to inclination effects, i.e., 
for pole-on sources v, » v sin /, and those are expected to be 
located mainly in the region v sin / < 150 km s ' . 

The stars 51 Oph, VV Ser, VX Cas, and SV Cep show 
<Wi()(Ha)> measurements that are lower than their correspond- 
ing <Wio(Ha)>,„,„ estimates. The reason could be uncertainties 
in the stellar parameters, introduced e.g. by close (unknown) 
companions or simply by the absence of Keplerian-rotating gas 
between the stellar surface and the co-rotation radius. In this case 
MA would hardly operate on these objects. Interestingly, all of 
them are again HBe stars according to their effective temper- 
ature (> 10000 K). Besides, the B-type object 51 Oph clearly 
deviates from the general trend in Fig. [8l likely pointing to an 
inner disk gas decoupled from the central star The nature of 
this object is not clear, since it could be a more evolved star or 
a y6-Pictoris analogue surrounded by a gas-rich debris disk (see 
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Fig. 9. (Top): Mean Ha profile of BF Ori from the October cam- 
paign in Paper I (solid line) and magnetospheric accretion syn- 
thetic profile that best reproduces it (dashed line; M, - 2.5 M0, 
R, = 2.5 R0, T, = 9000 K, magnetosphere size 4.2-5R,, M^cc 
= 10"*^ Mq yr"', i - 75 deg, and vsin/ - 37 km s"'). The dot- 
ted line is the same synthetic profile with no rotation. (Bottom): 
magnetospheric accretion synthetic Ha profiles with a stellar ro- 
tation rate of 75 km s"' (dashed lines) and with no rotation (dot- 
ted lines). 

e.g. lvan den Ancker et a lJ2001l:[StariretalJ2009l and references 
therein). 

In summary, we argue that the high rotational velocities of 
HAe stars mean that <Wio(Hq')> cannot be used as tracer of ac- 
cretion for those objects. At the same time, that the Ho- line width 
reflects the stellar rotation can be interpreted as a consequence 
of MA operating in the late-type HAeBe sources. 

5. Variability of accretion rates and accretion 
tracers 

The multi-epoch photometry and simultaneous spectra on which 
this work is based lead us to analyse the accretion rate variability 
and its possible relation with the changes in the accretion tracers. 
The photometry and spectra were taken during four campaigns 
of se veral days each, on different months (see Oudmaiier et al. 
12001 , and Paper I). The multi-epoch values for the Balmer ex- 
cess were derived from the individual U - B and B - V data 
taken on the same nights, using an expression analogous to that 
for <ADb> in Sect.O The multi-epoch Ha and [Oi]6300 lu- 
minosites were derived in paper 1 from line equivalent widths 
and the photometric measurements in the R band taken on the 
same nights. The line luminosities were dereddened using si- 
multaneous values for B - V and a procedure analogous to the 
one explained in Sect. |2]for the mean line luminosities. 

Most stars in our sample tend to show a constant Balmer 
excess, within the uncertainties. The strongest variations are be- 
tween 0.1 and 0.2 magnitudes, which are detected in some ob- 



servations of ~ 41% of the stars. An upper limit for the accre- 
tion rate variability can be estimated from the stars showing 
the smallest mean Balmer excesses and the largest Balmer ex- 
cess variations (see the right panel of Fig.[TJ. As extreme exam- 
ples, VI 686 Cyg shows an increase in the Balmer excess from 
0.04 magnitudes on 1998 May 15 to 0.18 magnitudes on 1998 
October 27. That change implies an accretion rate increase of 
a factor < 5. The Balmer excess of WW Vul decreased from 
0.14 magnitudes on 1998 May 15 to 0.04 magnitudes on 1998 
October 24, which translates into an accretion rate decrease of 
a factor < 4. Considering all the stars, our model yields a typ- 
ical upper limit of ~ 0.5 dex for the accretion rate changes, on 
timescales of days to months. More accurate constraints on pos- 
sible accretion rate variations below our detection limit would 
require UV spectra or narrow-band photometry. 

The analysis of the simultaneous variability of the Balmer 
excesses and the line luminosities is not only limited by the error 
bars provided for our broad-band photometry and mid-resolution 
spectra, but also by the number of data sampled per object. There 
is a maximum of 7 simultaneous photometric and spectral points 
per star, being typically 3 (see Table 3 in Paper 1). However, we 
find objects showing changes in the Balmer excess for which the 
corresponding Ho' and [Oi]6300 luminosities do not vary ac- 
cordingly, as would be expected for an accretion tracer. As an 
example, V1686 Cyg shows a Balmer excess close to mag- 
nitudes on 1998 May 15, which increased to the typical value 
from our averaged data (0.12 magnitudes) the following night. 
However, the corresponding Ha and [O i]6300 luminosities did 
not change. In contrast, there are stars showing variations in the 
line luminosities with a constant value of the Balmer excess. The 
one shown by BH Cep was ~ magnitudes during all observa- 
tions, but the Ha and [O i]6300 luminosities changed by a factor 
of almost 2. Two other stars showing a decoupled behaviour be- 
tween the Balmer excess and the line luminosities are provided 
in Fig. [To] RR Tau increased the Balmer excess from ~ 0.32 
magnitudes on 1998 October 27, to 0.45 magnitudes on 1999, 
January 30. Simultaneously, the Ha luminosity decreased by a 
factor 1 .4 and the [O i]6300 luminosity remained constant within 
the uncertainties -note that the slight variation in the equivalent 
width was not accompanied by a change in the line luminosity 
(Paper I)-. HK Ori showed a constant Balmer excess of ~ 0.57 
magnitudes at the same time that the Ha and [O i]6300 luminosi- 
ties decreased by a factor ~ 2 from 1998 October 25 to the fol- 
lowing night. The examples in this sense are not scarce, but the 
inverse case is not true; we do not detect stars in our sample for 
which a significant increase or decrease in the excess is reflected 
by the corresponding one in the line luminosities, but again, we 
note that most possible variations fall below our detection limit. 

Even considering the limitations of our data to analyse such 
small changes, the results described above suggest that the Ha 
and [Oi]6300 luminosities do not trace the accretion rate vari- 
ability in several stars, unless they reflect the accretion rate 
changes with a certain lag not covered by our data. Night- 
to-night variations in the atmospheric transmission could con- 
taminate the estimates of the variability in the Balmer excess. 
However, a careful photometric calibration was made, especially 
in the U band, by observing standard stars with different spectral 
types and airmasses and mak ing several tests after the calibra- 
tions (lOudmaijer et al.ll200ll) . It is unknown whether the varia- 
tions in other empirical tracers not analysed in this work could 
coiTelate with the accretion rate variability from the changes in 
the Balmer excess, or whether the variability could be aff'ected 
by other factors such as stellar companions. These topics are the 
subject of ongoing investigation. Multi-epoch, high-resolution 
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Fig. 10. Ha (left panels) and [O i]6300 (right panels) normalized 
profiles for RR Tau (top) and HK Ori (bottom), with different 
scales depending on the panel. The dashed and solid lines repre- 
sent the first and second observing dates (see text). 

spectra covering simultaneously the wavelength range from the 
UV to the near-lR will be very valuable in this sense. Those can 
be taken with the X-Shooter spectrograph on the VLT. 

6. Discussion 

The origin of the empirical calibrations between accretion and 
line luminosities could be explained from the influence that ac- 
cretion has on the line formation: Bot h the Ha and Bry lines 
have been explained from MA models (Muzerolle et al.|[l998bi 
1200 iL MDCH04), although the influence of winds can be non- 
negUgible (Kuros awa et al.ll2006i: fKraus et al.l l2008). The latter 
contribution seems to be more important as the stellar mass 
increases (see Paper I and references therein), which would 
explain the possible decrease in the slope of the Ho- calibra- 
tion, when compared to that for lower mass stars (left panel 
in Fig. |6]l. Accretion -powered outflow's could be the origin of 
the [Oi]6300 feature (ICorcoran & Ravl[T997l[T99l . They could 
therefore link the luminosity of this hne with that of accretion 
(middle panel in Fig. |6]l. If the [Oi]6300 emission is caused 
by the UV flux illuminating the disk surface ( Acke et al. 2005), 
then the UV excess from the accretion shock would explain the 
trend with the [Oi]6300 luminosity. However, several HAe stars 
with strong Balmer excesses do not show [Oi]6300 emission, 
regardless of their different degrees of disk flaring according to 
the Meeus et al. (2001) groups (see e.g. the group II object HD 
150193 and flie group I star V346 Ori). 

Although the mean Ha and [Oi]6300 luminosities are re- 
lated to the "typical" accretion rates (Sect. |4]i, the numerous 
cases for which the changes in the line luminosities and in the 
simultaneous Balmer excesses are decoupled suggest that those 
lines do not to trace the accretion rate variability. This could be 
the case for several other lines because the typically low accre- 
tion rate changes we have derived -less than 0.5 dex- contrast 
with the complex variability behaviour shown by the spectral 
tracers. Multi-epoch observations indicate that their variability 
amplitudes are not correlated: the Hei5876 line shows equiva- 
lent width variations over one order of magnitude for a given 
HAe star, at the sam e time that Ha varies typically by less than 
a factor 4 (Paper I). iNguven et all 12009) reports that the Can 
fluxes provide typical accretion rate variations of 0.35 dex while 



the Ha 10% width suggests changes of 0.65 dex for the same 
sample of CTT stars. These differences could be extended to a 
wide number of empirical tracers of accretion covering a wide 
wavelength range, indicating that the origin and strength of sev- 
eral of these lines are most probably influenced by different pro- 
cesses apart from accretion. 

These results make it reasonable to doubt that several empir- 
ical calibrations between the "typical" accretion and line lumi- 
nosities could not be caused by the accretion contribution to the 
lines but by a different source. We tentatively suggest that a com- 
mon dependence of both the accretion and the line luminosities 
on the stellar one could have strong influence, thereby driving 
the empirical calibrations. First, there seems to be a direct de- 
pendence of Lace on L» (Fig.|3] right). Second, it seems from the 
comparison between the data for low-mass stars and our results 
(see e.g. Herczeg & Hillenbrand 2008, and Fig. ^ that several 
line luminosities are also correlated with L»; e.g., most low-mass 
stars show Ha luminosities in the range -6 < logiLna I Lq) < 
-2, orders of magnitude lower than those shown by our more lu- 
minous stars. In addition, the line luminosities of the objects for 
which the shock models are not able to reproduce the Balmer 
excesses (see Sect. I3.2l i tend to be the highest in our sample, 
as are their stellar luminosities (see Table [T]i. From a statisti- 
cal point of view, the Spearman's probability (see e.g. lConoverl 
1980) that our accretion and Ha luminosities are uncorrected in 
Fig.|6]is only p = 4.1x10"'*%. However, this value does not con- 
sider the possible common dependence o f both luminosities on 
L*. The partial correlations technique (e.g.l Wall & Jenkinsl2003h 
is useful for this task, since it measures the degree of correlation 
between two variables, nullifying the effect of a common depen- 
dence on a third one. Under the hypothesis that both Lace and 
L(Ha) depend on L», the Spearman's probability from partial 
correlations increases to 27%. This value is close to that obtained 
from the comparison of two sets of random variables. A similar 
conclusion is reached for the Bry calibration and, to a lesser ex- 
tent, for the [O i]6300 one. These results indicate that the Lace 
- L/,„(. relations become much less significant once the common 
dependence on L* is removed. In this view, the apparent decrease 
in the slopes in the Lace - ^iwe trends for the intermediate-mass 
regime (Fig. |6ll would be reflecting the possible decrease in the 
Lace - L, correlation (right panel in Fig.|3]l. 

7. Conclusions 

We applied shock modelling within the context of MA to 38 
HAeBe stars, reproducing the strength of the observed Balmer 
excesses and deriving accretion rate estimates. The typical mass 
accretion rate is 2 x 10"^ Mq yr"', with scatter and a steep de- 
pendence on the stellar mass. The latter can be explained by 
the most massive HAeBe stars being the youngest, hence the 
strongest accretors. We obtained empirical expressions to relate 
accretion and the Ha, [O i]6300, and Bry luminosities, which are 
similar to those for the lower mass regime. However, there could 
be a slight decrease in the slope of the trends relating the accre- 
tion luminosity with both the line and the stellar luminosities. In 
contrast, the Ha line width at 10% of peak intensity can usually 
not be used to estimate accretion rates for HAeBe stars, unlike 
in the lower mass regime. The Ha width of HAe objects broad- 
ens as the projected rotational velocities increase, which is also 
consistent with MA operating in these stars. 

Although the "typical" accretion rates show clear trends with 
the mean line luminosities, the accretion-rate changes from the 
Balmer excess -typically lower than 0.5 dex- seem to be uncor- 
rected with the variability of the Ha and [Oi]6300 lines. We 



12 



I. Mendigutfa et al.: Accretion rates and accretion tracers of Herbig Ae/Be stars 



tentatively suggested that the origin of the empirical calibrations 
between the accretion and line luminosities could not be driven 
by the influence of accretion on the emission lines, but by a com- 
mon dependence on the stellar luminosity. The statistical level of 
significance of the empirical calibrations is strongly diminished 
once this dependence is considered. 

The shock model used is not able to reproduce the strong 
Balmer excesses shown by the four hottest stars in our sam- 
ple under reasonable input parameters. In addition, the Ha 10% 
width of four HBe objects is not consistent with the pres- 
ence of Keplerian gas between the stellar surface and the co- 
rotation radius, which is a major requirement for MA to operate. 
These results give additional support for a change from mag- 
netically driven accretion i n HAe stars to some other kind of 
accretion in HBes (see e.g. IVink et al]|2002l: [Eisner et al.ll2004l: 
[Mottram et al. 200%. 

It is finally pointed out that, although our results were ob- 
tained and interpreted within the context of MA, this means nei- 
ther that the authors found a way to test magnetospheric accre- 
tion onto HAeBe stars nor that the observational data could not 
be explained from a different scenario. In this respect, further 
observational and theoretical eff'orts on the magnetic fields nec- 
essary to drive MA in these objects are needed to establish firmer 
conclusions. 
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